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SUMMARY 


The  indie ial  lift  and  moment  functions  are  determined  approximately 
for  sinking  and  pitching  motion  at  Mach  numbers  M  of  0.5  and  0.6. 

These  functions  are  determined  from  a  knowledge  of  the  existing  oscil-  ■ 
latory  coefficients  at  the  low  reduced  frequencies  and  from  approximate 
expressions  of  these  coefficients  at  the  high  reduced  frequencies. 

The  beginning  portion  of  the  indicial  lift  function  associated 
with  an  airfoil  penetrating  a  sharp-edge  gust  in  subsonic  flow  is 
evaluated  by  use  of  an  exact  method.  By  use  of  an  approximate  method 
for  determining  the  remaining  portion,  the  complete  indicial  gust  func¬ 
tion  is  determined  for  M  =  0.5,  M  =  0.6,  and  M  =  0.7. 

All  the  indicial  lift  and  moment  functions  are  approximated  by  an 
exponential  series;  the  coefficients  which  appear  in  the  exponential 
approximations  for  each  indicial  function  are  tabulated  for  M  =  0.5, 

M  =  0.6,  and  M  =  0.7. 


INTRODUCTION 


In  the  study  of  transient  flows,  two  types  of  airfoil  motions  have 
had  special  significance  -  a  harmonically  oscillating  airfoil  and  an 
airfoil  experiencing  a  sudden  change  in  angle  of  attack.  The  lift  and 
moment  functions  for  an  airfoil  experiencing  a  sudden  change  in  angle 
of  attack  are  commonly  referred  to  as  indicial  lift  and  moment  func¬ 
tions.  In  references  1  and  2  the  indicial  lift  and  moment  functions 
of  a  two-dimensional  airfoil  for  sinking  and  pitching  motion  in  a  sub¬ 
sonic  compressible  flow  are  presented  for  M  =  0.7.  ■  The  functions  are 
determined  from  the  reciprocal  relations  between  the  indicial  functions 
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and  the  oscillatory  coefficients.  Examination  of  the  oscillatory  coef¬ 
ficients  calculated  by  Dietze,  Schade,  and  Possio  (see,  for  example, 
refs.  3,  4,  and  5)  indicates  that  data  are  not  available  over  a  suffi¬ 
cient  range  of  reduced vfrequencies  for  the  direct  utilization  of  this 
method  to  determine  the  indie ial  lift  functions  at  other  subsonic  Mach 
numbers . 

The  indicial  functions  for  M  =  0.8  have  been  evaluated  by  Lomax, 
Heaslet,  and  Sluder,  in  reference  6,  employing  an  alternate  method 
applicable  for  both  subsonic  and  supersonic  Mach  numbers.  The  alter¬ 
nate  method  for  subsonic  Mach  numbers  is  characterized  by  an  exact, 
readily  evaluated  expression  for  the  beginning  portion  of  the  indicial 
functions  and  an  approximate  expression  for  the  remaining  portion.  The 
solution  for  the  remaining  portion  of  the  indicial  function  is  lengthy 
and  tedious,  and,  consequently,  an  alternate  approach  is  used  in  the 
present  paper  to  calculate  this  portion  of  the  indicial  functions. 

In  the  present  paper  indicial  functions  for  M  =  0.5  and  M  =  0.6 
are  obtained  in  approximate  form  for  an  extensive  range  of  airfoil 
travel.  The  beginning  portions  of  the  indicial  functions  are  calcu¬ 
lated  by  the  method  of  reference  6.  The  difficulties  of  calculating 
the  remaining  portions  of  the  indicial  functions  by  the  method  of 
reference  6  are  avoided  by  approximating  the  missing  portions  of  the 
oscillatory  coefficients  of  references  3  and  k  so  that  the  method  of 
references  1  and  2  could  be  used.  The  approximation  is  made  in  terms 
of  the  initial  parts  of  the  indicial  functions  and  the  known  parts  of 
the  oscillatory  coefficients.  An  estimate  of  the  reliability  of  the 
indicial  functions  is  also  presented. 

During  the  last  stage  of  preparing  the  paper,  it  was  found  that 
Timman,  Van  de  Vooren,  and  Greidanus  have  computed  the  oscillatory 
coefficients  in  terms  of  known  functions  for  Mach  numbers  0.35*  0.5, 

0.6,  0.7,  and  0.8  and  for  reduced  frequencies  ranging  from  0.1  to  3*0 
(see  ref.  7)«  A  discussion  of  the  effects  of  differences  between  the  « 
data  of  reference  7  and  that  of  the  present  investigation  on  the  results 
of  the  calculations  is  included. 


SYMBOLS 


s  distance  traveled,  half-chords 

a)  angular  frequency 

V  forward  velocity  of  airfoil 


c 


chord 


> 


•4 


* 
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k 

*0 

KkQ 

^00 


€max 


L(s) 

M(s) 


P 

h 

0 

<1 


k]_(s ) 


®l(a) 


aq(s) 


3 

Cartesian  coordinates 

reduced- frequency  parameter,  ox:/2V 

highest  reduced  frequency  for  which  oscillatory 
coefficients  are  known 

value  of  reduced  frequency  where  maximum  error  in 
approximated  oscillatory  coefficients  occurs 

error  in  oscillatory  coefficients  at  infinite  reduced 
frequency 

maximum  error  in  oscillatory  coefficients  for 
region  K  >  1 

section  lift,  positive  in  downward  direction 

section  pitching  moment  about  quarter-chord  point, 
positive  when  moment  tends  to  depress  trailing 
edge 

dens ity 

amplitude  of  vertical  or  sinking  displacement, 
positive  downward,  half-chords 

angle  of  pitch,  positive  when  trailing  edge  is 
lower  than  leading  edge 

rate  of  change  of  angle  of  pitch  with  respect  to 
distance  traveled  in  half-chords,  positive  when 
trailing  edge  is  falling  with  reference  to  leading 
edge 

indicial  lift  function  for  a  two-dimensional  air¬ 
foil  experiencing  a  sudden  change  in  vertical 
velocity 

indicial  moment  function  for  a  two-dimensional  air¬ 
foil  experiencing  a  sudden  change  in  vertical 
velocity  where  moment  is  taken  about  quarter- 
chord  point 

indicial  lift  function  for  a  two-dimensional  air¬ 
foil  experiencing  a  sudden  change  in  pitching 
velocity  about  its  leading  edge 


M 

CP 

mlq(s) 

U 

k2(s) 

Fc(k)  +  iGc(k) 
M(k)  +.  iN(k) 
FCq(k)  +  iGcq(k) 
Mq(k)  +  urq(k) 

f(k)  =  Fc(k)  -  Fc 
ra(k)  =  M(k)  -  M(oo 
fq(k)  =  FCq(k)  -  1 

mq(k)  =  Mq(k)  -  M( 

an>^n^an^^)n 

A,B 


Mach  number 

pressure  coefficient 

indie ial  moment  function  for  a  two-dimensional  air¬ 
foil  experiencing  a  sudden  change  in  pitching 
velocity  about  its  leading  edge,  moment  taken 
about  quarter-chord  point 

velocity  of  sharp-edge  gust 

indicial  lift  function  for  a  two-dimensional  air¬ 
foil  due  to  penetration  of  a  sharp-edge  gust 

complex  compressible  oscillatory  lift  derivative 
due  to  sinking  velocity 

complex  compressible  oscillatory  moment  derivative 
due  to  sinking  velocity 

complex  compressible  oscillatory  lift  derivative 
due  to  pitching  velocity  only 

complex  compressible  oscillatory  moment  derivative 
due  to  pitching  velocity  only 


oo) 


Cq(“) 

(oo) 


constants  in  exponential  approximations 
constants  defined  by  equations  (13)  and  (l4) 


Subscript: 

3c/4 


for  airfoil  rotating  about  three- quarter-chord  point 
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METHOD  OF  ANALYSIS 


The  indicial  functions  are  determined  in  two  parts.  In  the  region 
0  =  s  S  ■  ,  readily  determined  indicial  functions  given  in  refer¬ 
ence  6  are  used.  For  the  region  s  >  — ‘ ^  ,  the  relatively  more  dif- 

ficult  method  of  reference  6  is  avoided  by  an  alternate  approach  which 
utilizes  the  reciprocal  relations  between  the  indicial  and  the  oscil¬ 
latory  functions  given  in  references  1  and  2.  Before  these  relations 
could  be  applied  at  M  =0.5  and  M  =  0.6,  it  was  necessary  to  deter¬ 
mine  the  unknown  oscillatory  coefficients  at  the  high  reduced  frequen¬ 
cies,  k  >  1.  Approximate  expressions  for  these  coefficients  are 
obtainable  by  combining  the  existing  data  for  the  oscillatory  coeffi¬ 
cients  at  the  low  reduced  frequencies,  0  <  k  ^  1,  with  the  indicial 

oM 

functions  for  the  region  0^3^  — -  through  the  use  of  reciprocal 

1  +  M 

relations.  The  reciprocal  relations  are  then  used  to  obtain  the 
remaining  portion  of  the  indicial  functions. 

The  following  procedure  was  used  to  determine  the  oscillatory  coef¬ 
ficients  at  the  high  reduced  frequencies:  The  reciprocal  relation 
between  the  indicial  lift  function  due  to  sinking  velocity  and  the  in- 
phase  oscillatory  lift  coefficient  is  (see  eq.  (lh)  of  ref.  1) 


*!(»)  =  FCM  *  I  r  f(t>kln  “  ■  <»  >  0)  (la) 

or  alternate ly, 

r  fOOgjgJg.  ik  .  |  ti(s)  .  g  T  (.)  .  °  1-Mam  ka  dt  (lb) 

ko  Jo  k 

where  f(k)  =  Fc(k)  -  Fc(°°). 

In  equation  (lb)  values  of  k^(s)  are  known  from  the  method  of 
reference  6  for  0  ^  s  ^  The  value  of  Fc(«)  which  has  been 

shown  in  reference  1  to  correspond  to  the  value  of  k]_(s )  at  s  =  0 
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is }  consequently,  also  known.  The  oscillatory  coefficient  f(k)  is 
considered  known  for  the  range  of  reduced  frequencies  0  <  k  <  kQ. 
Thus,  the  right-hand  side  of  equation  (lb)  is  known  for  values  of 


0  <:  s  ^ 


2M 

1  ’+  M* 


In  the  region  kQ  <  k  <  00,  f(k)  is  unknown  and  is 


to  be  determined. 


A  solution  for  f(k)  in  the  region  kQ  <  k  <  00  can  be  expressed 
by  a  finite  convergent  series  with  unknown  coefficients.  Definition 
of  these  coefficients  is  possible  by  a  simultaneous  solution  of  a  set 
of  algebraic  equations.  The  set  of  equations  is  formulated  by  sub¬ 
stituting  the  finite  series  for  f(k)  (together  with  appropriate 
boundary  conditions)  in  the  integral  on  the  left  side  of  equation  (lb). 
This  integral  is  evaluated  to  obtain  an  algebraic  expression  in  terms 
of  the  unknown  series  coefficients  and  the  variable  s.  Each  simul¬ 
taneous  equation  in  the  set  is  then  obtained  by  assigning  to  the 
variable  s  a  value  selected  in  the  range  where  k]_(s)  is  known. 

The  number  of  these  equations  should  correspond  to  the  number  ot 
series  coefficients  to  be  determined.  The  set  of  equations  is  then 
solved  simultaneously  for  the  series  coefficients  in  the  approximation 
for  f(k). 


The  accuracy  of  the  approximate  solution  for  f(k)  depends  on 
the  range  of  s  over  which  k]_(s)  is  known,  the  value  of  kQ,  and 
the  number  of  terms  in  the  series.  The  error  introduced  by  a  finite 
series  for  f(k)  is  minimized  by  the  fact  that  the  oscillatory  coef¬ 
ficient  for  the  higher  reduced  frequencies  is  associated  primarily 
with  the  initial  portion  of  the  indicial  function  which  is  known. 


Once  the  approximate  expression  for  f(k)  is  determined  at  the 
high  reduced  frequencies,  k  >  1,  this  expression  and  the  oscillatory 
coefficients  at  the  low  reduced  frequencies,  0  <?  k  ^  1,  can  be  sub¬ 
stituted  into  equation  (la)  to  determine  kq(s)  for  values  of 
2M 

s  >  - . 

1  +  M 


The  same  over -all  method  may  be  used  for  determining  the  indicial 
moment  on  a  sinking  airfoil  as  well  as  the  indicial  lift  and  moment 
functions  for  a  pitching  airfoil.  The  reciprocal  relation  used  for 
approximating  the  oscillatory  moment  coefficients  at  the  high  reduced 
frequencies  due  to  sinking  motion  is  as  follows: 


* 


* 
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where  m(k)  is  defined 

m(k)  =  M(k)  -  M(oo) 


For  the  case  of  a  pitching  airfoil,  it  was  shown  in  reference  2  that 
the  only  additional  indicial  function  that  need  he  considered  is  the 
one  associated  with  the  pitching  velocity  of  the  airfoil.  The  recip¬ 
rocal  relations  used  in  this  analysis  for  detennining  the  oscillatory 
lift  and  moment  coefficients  at  the  high  reduced  frequencies  as  well  as 
the  complete  indicial  functions  for  the  case  of  pitching  velocity  are 
for  the  lift,  ’ 


klq(s)  =  FCq(“) 


fq(k)sin  ks 
k 


dk 


(s  >  0)  (3) 


where 

fq  =  Fcq00  -  %U)  (4) 

and,  for  the  moment. 


mq(k)sin  ks 
k 


dk 


(s  >  0)  (5) 


where 


mq(k)  =  Mq(k)  -  Mq(oo) 


(6) 


For  convenience,  a  summary  of  the  reciprocal  relations  for  lift  and 
moment  due  to  both  sinking  and  pitching  motion  is  given  in  table  I. 
together  with  the  definition  of  the  functions  appearing  in  the 
reciprocal  relations. 
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AVAILABLE  OSCILLATORY  COEFFICIENTS  AND  THEIR  APPROXIMATION 
AT  HIGH  REDUCED  FREQUENCIES 


Summary  of  available  oscillatory  coefficients  at  M  =  0.5  and 
M  =0.6.-  The  two-dimensional  compressible  oscillatory  coefficients 
for  the  real  and  imaginary  parts  of  the  lift  and  moment  for  sinking 
and  pitching  motion  are  given  in  table  II  for  M  =  0.5  and  M  =  0.6. 
The  results  given  in  table  II  were  taken  from  two  sources  (see  refs.  3 
and  4)  and  have  been  converted  to  the  form  given  by  the  equations 
noted  in  the  last  column  of  table  I  for  the  lift  and  moment.  The 
range  of  reduced  frequencies  taken  from  each  source  is  indicated  in 
table  II.  Expressions  for  the  end  points  of  the  oscillatory  coeffi¬ 
cients,  that  is,  for  k  =  0  and  k  =  00,  were  determined  from  the 
expressions  given  in  references  1  and  2. 

Expressions  for  indie ial  lift  and  moment  functions  in  region’ 

0  <  a  ^  — .  -  The  expressions  for  the  indicial  lift  and  moment  func- 

_ 1  +  M 

tions  due  to  sinking  and  pitching  motion  that  are  given  in  reference  6 

for  the  region  0  ^  s  'S  -  have  been  converted  to  the  notation 

1  +  M 

used  in  this  paper  and  are  as  follows: 


ki(s) 


(7) 


m 


l(s)  = 


2 

IrtK 


M(1  ■  M)  ■  W2  ■  M) 


(8) 


klq(s) 


1 

jtM 


1  "  2M^  "  M)  +  4mV 


Sid  .  M\ 

liMV  ?  / 


(9) 


mp  (s)  =  -  — 

V  JtM 


-  1  £(l 
8  MU 


M) 


32M‘ 


(2  -  6M  +  3M2)  + 


12M3  L; 


i(l  -  M)3  + 


(10) 


MS 
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Approximation  of  oscillatory  coefficients  at  high  reduced  fre¬ 
quencies.-  The'  following  series  was  assumed  as  an  approximate  expre, 
sion  for  f(k)  in  the  range  kQ  ^  k  ^ 


f(k)  = 


-H 


ke-an(k-ko) 


where  an  and  are  the  unknown  coefficients  to  be  determined  for 

f(k)  at  a  given  Mach  number.  Although  several  forms  of  series  could 
be  used  to  express  f(k),  this  form  was  found  to  be  the  most .convenient 
for  purposes  of  integration.  Note  that  equation  (ll)  satisfies  the 
condition  that 


lim  f(k)  =  lim 
k-*  co  k  — ^  00 


|jc(k)  -  Fc(»j]  = 


Two  boundary  conditions  which  require  the  approximate  solution 
for  f(k)  in  the  region  kQ  =  k  5-  °°  to  join  the  known  value  at  k 
continuously  in  magnitude  and  slope  are 


f(k) 
.  k 


These  two  conditions  are  used  to  reduce  the  number  of  unknowns  in 
equation  (ll).  Substitution  of  these  two  conditions  into  equation  (ll) 
leads  to  the  following  expressions  for  the  coefficients  a^  and  a2: 

B  +  Aap 

ai  - - -  (15) 

cep  -  a]_ 


a]_A  +  B 
al  "  a2 


(16) 
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If  the  expression  for  f(k)  given  by  equation  (ll)  is  substituted 
into  the  left-hand  side  of  equation  (lb)  and  if  the  values  of  aq  and 
ag  are  replaced  by  the  expressions  given  by  equations  (15 )  and  (l6), 
then  two  expressions  from  which  aq  and  a2  can  be  evaluated  are 
obtained.  For  the  case  M  =  0.5  and  M  =  0.6,  the  value  of  k q  is 
unity,  and  the  two  desired  expressions  are 


B  +  Aa2  Sq  cos  sq  +  aq  sin  Sq 

_____  —  -  + 


aqA  +  B  s^  cos  Sq  +  a2  sin  Sq 
al  "  a2  +  sq^ 


f(k)  sin  ksq 
k 


dk 


(17a) 


and 


B  +  Aa2  S2  cos  S2  +  aq  sin  S2  aqA  +  B  Sp  cos  s2  +  a2  sin  s2 
a2  ~  al  aq2  +  s22  al  "  a2  ^2  +  s22 

!&<■*>  -  (17b) 

where  ai  and  aP  are  numbers  greater  than  zero,  and  s-i  and  sP 

2M* 

are  greater  than  zero  but  less  than  or  equal  to  - . 

1  +  M 

A  simultaneous  solution  of  the  nonlinear  equations  (17a)  and  (17b) 
for  the  coefficients  aq  and  a2  was  obtained  by  a  graphical  procedure. 

In  some  cases  only  an  approximate  solution  for  aq  and  a2  could  be 
obtained.  In  these  instances,  values  for  aq  and  a2  were  selected 
such  that  the  magnitudes  of  error  incurred  in  the  indicial  function 
were  equal  at  the  points  sq  and  s2.  For  the  case  of  M  =  0.5,  values 

of  sq  =  O.lOn  and  S2  =  0.20jt  were  chosen  and,  for  M  =  0.6, 
sq  =  0.12n  and  s2  =  0.23*t  were  chosen.  Both  of  these  values  are 

2M 

within  the  limits  0  <  s  <  -  and  are  approximately  equally  spaced. 

1  +  M 
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Values  for  the  constants  ax  and  a2  could  then  he  determined  from 
equations  (15)  and  (l6).  In  table  III  values  for  the  constants 
a2)  al>  and  <x2  which  appear  in  the  exponential  approximation 

(eq.  (11))  are  summarized  for  the  lift  and  moment  due  to  both  sinking 
and  pitching  motion. 

RESULTS  AND  DISCUSSION 


Determination  of  indicial  functions.-  With  the  coefficients  in 
equation  (11)  determined  for  the  lift  and  moment  oscillatory  coeffi¬ 
cients  due  to  both  sinking  and  pitching  motion,  the  indicial  functions 
kp(s),  m^(s),  k-^  (s),  and  m-^  (s)  were  evaluated  over  the  range 

2M  ^  Q  . 

S  >  T  +  M  for  M  =  °-5  and  M  =0.6  by  the  methods  of  references  1 
and  2.  These  results  and  corresponding  functions  calculated  by  the 

method  of  reference  6  for  Oi  s  ^  ^  are  shown,  respectively,  in 

figures  1,  2,  3,  and  k  together  with  the  results  for  M  =  0.7,  obtained 
from  references  1  and  2.  Also  included  is  the  function  for  incompressi¬ 
ble  flow  obtained  by  Wagner.  The  indicial  lift  and  moment  functions  due 
to  pitching  velocity  for  the  airfoil  rotating  about  the  thre e -quarts r- 
chord  point  are  shown  in  figures  5  and  6.  For  incompressible  flow, 
the  lift  and  moment  are  impulsive  at  s  =  0  and  constant  thereafter; 
however,  for  compressible  flow,  the  lift  and  moment  are  finite  and 
time -dependent  and  approach  a  constant  value  very  rapidly.  Comparison 
of  the  indicial  lift  functions  due  to  both  sinking  and  pitching 
velocities  for  subsonic  Mach  numbers,  figures  1  and  3,  indicates  that 
as  the  subsonic  Mach  number  increases  the  lift  at  s  =  0  becomes  less 
and,  for  s  equal  to  more  than  approximately  two  half-chords,  the 
growth  of  lift  to  the  steady  state  becomes  less  rapid.  For  the  case 
of  the  moments  about  the  quarter- chord  point  due  to  sinking  and 
pitching  motion,  it  is  similarly  seen  from  figures  2  and  4  that  while 
the  moments  are  impulsive  at  zero  for  incompressible  flow  the  moments 
are  finite  and  their  values  decrease  with  increasing  Mach  number  for 
compressible  flow.  Although  the  moments  for'  s  >  0  are  constant  for 
incompressible  flow,  a  time -dependent  moment  which  rapidly  approaches 
a  constant  exists  for  compressible  flow  where  the  rapidity  of  approach 
decreases  with  increasing  Mach  number. 

Considerations  for  indicial  lift  due  to  penetration  of  a  sharp-edge 
gust.-  In  references  1  and  6  it  was  shown  that  the  indicial  lift  func- 

tion  k]_(s)  for  the  region  s  > -  is  associated  primarily  with 

1  -  M  . 

the  lift  due  to  circulation  only.  As  a  consequence,  an  approximation 
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for  the  indicial  lift  function  k2(s)  due  to  penetration  of  a  sharp- 
edge  gust  in  this  region  can  he  determined  on  the  basis  of  the  relation 
ship  that  exists  between  the  two  indicial  lift  functions  ki(s)  and 
k2(s)  for  the  circulatory  portion  of  the  lift  for  incompressible  flow. 

In  reference  1  the  function  k2(s)  is  determined  for  the  region 
s  >  _  on  the  basis  of  this  relationship  for  M  =  0.7.  However, 

no  information  concerning  the  beginning  portion  of  this  function  is 
available  at  present. 

A  method  for  determining  an  explicit  expression  for  this  function 

in  the  region  0  <  s  <  is  possible  from  an  extension  of  the 

analysis  given  in  reference  6  where  the  boundary  conditions  for  the 
perturbation  velocity  due  to  the  penetration  of  a  gust  can  be  sub- 
stituted  for  the  perturbation  velocity  due  to  sinking  motion.  Details 
of  the  evaluation  of  the  function  k2(s)  for  this  region  are  given  in 

appendix  A.  The  expression  obtained  for  this  function  is  relatively 
simple  as  can  be  seen  from  the  following  equation: 


k2(s) '  -7M 


(° 


(18) 


where  k2(s)  is  defined  such  that  the  lift  on  the  airfoil  is 

L(s)  =  -JtpcVUk2(s)  (19) 


The  gust  function  k2(s)  for  the  region  =  s  -  m  Can 

be  determined  from  an  extension  of  the  results  given  in  reference  6. 
An  explicit  expression  cannot  be  easily  obtained  because  of  the  com¬ 
plexity  of  the  downwash  regions  that  must  be  considered.  Since  the 
penetration  function  is  known  everywhere  except  in  the  region 

2**  <  s  <  — _  a  reasonable  estimate  of  the  curve  may  be  obtained 

2M 


“1-M 

by  joining  the  ends  of  the  two  solutions  at  the  points 


1  +  M 


and 


2M 


1  -  M 


by  a  smooth  curve. 


In  figure  7,  the  indicial  lift  functions  due  to  penetration  of  a 
sharp-edge  gust  are  shown  for  M  =  0.5>  M  =  0.6,  and  M  =  0.7  ^or  an 
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extensive  range  of  airfoil  travel.  Also  shown  is  the  solution  for  M  =  0 
obtained  by  Kussner.  Inspection  of  figure  7  indicates  that  not  only  is 
the  growth  of  lift  to  the  steady  state  less  rapid  as  the  subsonic  Mach 
number  is  increased  but  also  the  magnitude  of  lift  in  the  beginning 
portion  is  less. 

Approximation  of  indicia!  functions  by  analytic  expressions . -Since 
the  exponential  function  has  a  simple  operational  equivalent  and  a  series 
-of  such  functions  has  been  found  convenient  for  approximating  kp(s)  at 
M  =  0  and  M  =  0.7,  a  finite  sum  of  exponential  terms  was  chosen  to 
approximate  the  indie ial  lift  and  moment  functions  for  sinking  and  « 
pitching  motion  at  M  =  0.  5  and  M  =  0.6  as  follows: 

bo  +  bj_e~^ls  +  b2^  +  b^e  ^3S  (20) 


In  table  IV  values  for  these  constants  are  given  for  the  lift  and 
moment  due  to  sinking  motion  as  well  as  for  the  lift  and  moment  due  to 
pitching  motion  about  the  three-quarter-chord  point.  The  constants 
are  noted  for  M  =  0.5  and  M  =  0. 6  together  with  those  for  M  =  0.7 
from  reference  2.  In  reference  8  the  constants  are  given  for  incom¬ 
pressible  flow.  Plots  of  these  exponential  approximations  are  shown 
in  figures  1,  2,  5>  and  6.  Comparison  of  the  approximations  with  the 
actual  curves  in  these  figures  indicates  good  agreement. 

The  exponential  approximation,  given  by  equation  (20)  was  also 
found  to  be  convenient  for  fitting  the  indicial  lift  functions  k2(s) 
due  to  penetration  of  a  sharp-edge  gust  shown  in  figure  7.  In  table  IV 
the  constants  required  for  fitting  the  k2(s)  function  at  M  =  0.5, 

M  =  0.6,  and  M  =  0.7  tor  use  in  equation  (20)  are  given.  In  refer¬ 
ence  8  the  constants  for  this  indicial  function  are  given  for  incom¬ 
pressible  flow.  Comparison  of  these  approximations  with  the  derived 
curves  (see  fig.  7)  indicates  excellent  agreement. 

Approximation  of  oscillatory  coefficients  by  analytic  expressions. - 
The  corresponding  approximate  expressions  for  the  harmonically  oscil¬ 
lating  airfoil  can  be  found  from  equation  (20)  and  !an  alternate  form  of 
the  reciprocal  equations.  For  the  case  of  the  lift  due  to  sinking 
motion,  the  reciprocal  equation  is 


Fc(k)  +  iGc(k)  =  ik 


k]_(s ) 


-iks 


0 


ds 


(21) 
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Substitution  of  the  approximation  for  k]_(s)  given  by  equation  (20) 
leads  to  the  following  expressions  for  Fc(k)  and  Gc(k) 


Fc(k)  =  b0  + 


bj_k^  b2k2  b3k2 

Pi2  +  k2  p22  +  k2  po2  +  k2 


.  ,  blPlk  b2P2k 

Gc(k)  —  — - -  +  -  4- 


b3P3k 


Px2  +  k2  p22  +  k2  P32  +  k2 


In  a  similar  fashion,  expressions  for  the  moment  functions  M(k)  and 
N(k)  can  be  obtained  in  the  form  given  by  equation  (22)  as  well  as 
for  the  lift  functions  due  to  pitching  velocity  about  the  three- 
quarter-chord  point  (Fcq)3c/4(k)  and  (GCq)3c/4(k)  and  the  moment 

functions  (Mq)3C//i,.(k)  and  (Nq)3c/lj.(k) .  The  transfer  of  these  coef¬ 
ficients  from  the  leading  edge  to  the  three-quarter-chord  point  was 
made  with  the  aid  of  equations  (27)  of  reference  2. 

Estimate  of  reliability  of  indicial  functions.-  Although  the 

indicial  functions  are  known  exactly  in  the  region  0  <  s  <  —  the 

=  1  +  M' 

reliability  of  these  functions  for  values  of  s  >  — — —  is  question- 

able  since  the  oscillatory  coefficients  at  the  high  reduced  frequencies, 
which  were  used  in  determining  the  indicial  function  in  this  region 
have  been  determined  only  approximately.  An  estimate  of  the  reliability 
of  the  indicial  functions  in  this  region  due  to  the  errors  in  the  oscil¬ 
latory  coefficients  at  the  high  reduced  frequencies  could  be  made  on 
the  basis  of  the  differences  which  exist  between  the  approximate  solu¬ 
tion  (based  on  the  oscillatory  coefficients)  and  the  exact  solution 

(based  on  ref.  6)  of  the  indicial  functions  in  the  region  0  <  s  <  — 

o-  ~  1  +  M' 

Since  the  oscillatory  coefficients  at  the  high  reduced  frequencies 

principally  affect  the  beginning  portion  of  the  function  kq(s),  the 
differences  between  the  approximate  and  exact  solutions  in  this  region 
should  give  a  reasonable  and  conservative  estimate  of  the  errors 
incurred  in  the  indicial  functions  for  larger  values  of  s. 


MCA  TN  2739 


15 


In  appendix  B,  a  chart  (see  fig.  8)  is  obtained  for  indicating  the 
v&Fi&t ion  of  errors  in  the  indicial  functions  with  airfoil  travel  s 
due  to  the  errors  in  the  in-phase  oscillatory  coefficients  at  the  high 
reduced  frequencies.  Although  figure  8  indicates  the  errors  in  only 
the  indicial  lift  function  k^(s),  the  results  in  this  figure  are  also 
applicable  for  the  other  indicial  functions.  Inspection  of  figure  8 
for  kQ  =  1  indicates  that  the  maximum  error  in  the  indicial  function 
occurs  within  a  quarter  of  a  chord  regardless  of  the  position  of  the 
maximum  error  occurring  in  the  in-phase  oscillatory  coefficients  (as 
indicated  by  several  values  of  k).  Since  s  =  0. 5  lies  within  the 

region  0  <  s  <  +~ M  ^°r  ^  ~  ®.5  and  M  =  0.6,  it  may  be  expected 

that  the  greatest  errors  in  the  indicial  functions  occur  in  this  region 
and  that  errors  beyond  this  point  approach  zero  very  rapidly.  The 
largest  differences  between  the  indicial  functions  determined  from  the 
oscillatory  coefficients  and  by  the  method  given  in  reference  6  in  the 
2m 

region  0  <  s  <  - - -  were  found  to  be  less  than  1. 5  percent  for 

1  +  M.  ^ 

M  =  0.5  and  M  =  0.6.  Thus,  it  appears  from  figure  8  that  the  largest 

errors  in  the  indicial  functions  beyond  the  point  s  =  -  are  much 

1  +  M 

less  than  1. 5  percent  and  appear  to  approach  zero  very  rapidly  for 
increasing  values  of  airfoil  travel  s. 

Effects  of  discrepancies  between  data  of  reference  7  and  those  of 
present  investigation  on  indicial  functions. -  Just  prior  to  completion 
of  this  paper  oscillatory  lift  and  moment  coefficients  which  cover  a 
larger  range  of  reduced  frequencies  than  that  considered  by  Dietze  and 
Schade  were  published  (see  ref.  7).  In  figures  9  and  10  a  comparison 
is  shown  of  the  in-phase  coefficients  used  in  the  present  analysis 
with  those  obtained  from  reference  7*  For  the  range  of  reduced  fre¬ 
quencies  0  <  k  ^  2.5>  the  data  agree  within  approximately  3  percent. 

However,  for  k  >  2. 5  the  differences  tend  to  become  larger  than 
3  percent  in  some  cases. 

Consideration  of  the  over-all  effect  of  the  discrepancies  shown 
in  figures  9  and  10  indicates  that  these  differences  would  result  in 
a  maximum  error  of  3  percent,  in  the  indicial  functions.  This  amount 
of  error  does  not  appear  serious,  and  in  view  of  the  few  data  avail¬ 
able  from  reference  7,  recomputation  of  the  indicial  function  does 
not  appear  worthwhile.  It  might  be  noted  that  the  data  of  Dietze  and 
Schade  agree  better  mutually  than  with  the  data  of  reference  7.  The 
reason  for  these  discrepancies  is  not  known;  consequently,  the  selection 
of  the  most  accurate  oscillatory  coefficients  is  still  questionable. 
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CONCLUDING  REMARKS 


From  approximate  expressions  for  the  in-phase  oscillatory  coeffi¬ 
cients  determined  for  the  high  reduced  frequencies  and  with  a  knowledge 
of  the  existing  data  at  the  low  reduced  frequencies,  the  indie ial  lift 
and  moment  functions  have  been  determined  for  sinking  and  pitching 
motion  at  Mach  numbers  M  of  0.5  and  0.6. 

The  indie ial  lift  function  due  to  penetration  of  a  sharp-edge  gust 
is  determined  for  M  =  0.5,  M  =  0.6,  and  M  =  0.7  for  an  extensive 
range  of  airfoil  travel.  Inspection  of  these  functions  indicates  that 
not  only  is  the  growth  of  lift  to  the  steady  state  less  rapid  as  the 
subsonic  Mach  number  is  increased  but  also  the  magnitude  of  lift  in 
the  beginning  portion  is  less. 


Langley  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Langley  Field,  Va. ,  April  15,  1952. 
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APPENDIX  A 


DETERMINATION  OF  INDICIAL  LIFT  FUNCTION  k2(s)  DUE  TO  PENETRATION 

OF  A  SHARP -EDGE  GUST  FOR  0  <  s  <  , m 

"  1  +  M. 


In  reference  6  a  method  is  outlined  for  determining  the  subsonic 
indie ial  lift  functions  for  a  sinking'  airfoil.  However,  if  the  boundary 
conditions  for  the  perturbation  velocity  due  to  the  penetration  of  a 
gust  are  substituted  for  the  perturbation  velocity  due  to  sinking  motion, 
the  function  k2(s)  can  be  obtained.  Because  of  the  complexity  of  the 
problem  arising  after  the  Mach  forecone  intersects  the  trailing  edge, 

the  analysis  is  made  only  for  the  region  0  ^  s  %  - - The  boundary 

condition  for  the  gust  velocity  U  differs  from  the  boundary  condition 
of  the  sinking  airfoil  in  that  there  is  no  perturbation  velocity  on  the 
part  of  the  airfoil  left  unshaded  in  the  following  sketch: 
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If  as  in  reference  6  the  axes  in  the  preceding  sketch  are  transformed 
hy  the  relations 


E'if(S+x)  ;  n-^(s-x) 

then  the  geometry  can  he  represented  hy 


The  expression  for  the  pressure  coefficient  Cp  can  he  written  in 
terms  of  the  coordinate  system  as  follows,  where  the  subscripts 

indicate  the  region  for  which  the  particular  equation  applies: 

CPl  =  0  (Al) 


Here  it  might  he  noted  that  equation  (A2)  differs  from  the  corresponding 
equation  in  reference  6  only  hy  the  lower  limit  of  the  last  integral 
since  this  integration  extends  over  only  that  part  of  the  wing  which 
has  penetrated  the  gust. 
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APPENDIX  B 

INVESTIGATION  OF  ERRORS  IN  INDICIAL  FUNCTION  DUE  TO  ERRORS  IN  IN-PHASE 
OSCILLATORY  COEFFICIENTS  AT  HIGH  REDUCED  FREQUENCIES 

Consider  the  following  reciprocal  relation  for  the  lift  due  to 
sinking  velocity: 

(s)  _  2  f"  Fc(k)sin  to  dk  (B1) 

rt  Jq  k 

It  is  assumed  that  Fc(k)  is  known  for  values  of  0  ^  k  ^  ko;  however, 
for  values  of  k  >  kQ  an  error  in  the  Fc(k)  function  is  assumed  of 
the  following  form: 

Fc00error  =  ±6°o(l  "  t)  (B2) 

where  e^,  is  the  magnitude  of  the  error  at  k  = 
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If  equation  (B2)  is  substituted  into  equation  (Bl),  an  expression  for 
the  error  in  the  kp(s)  function  may  be  obtained  as  follows: 


=  ^  |  -  sin  kos  -  Si(kQs)  +  kQS  Ci(kQs)  (B3) 


where  Si(x)  and  Ci(x)  are  tabulated  in  Jahnke  and  Emde  (ref.  9)* 

If  the  value  of  Fc(°o)  is  known  and  if  the  maximum  error  is 
assumed  to  occur  at  some  arbitrary  point  rckQ  (where  k  >  1) ,  two 
error  functions  of  the  type  given  by  equation  (B2),  where  one  is  dis¬ 
placed  from  the  point  k q  by  an  amount  KkQ,  may  be  superposed  such 

that  the  value  of  Fc(co)error  =  0  as  illustrated  in  the  following 
sketch: 

Pc(k) 


The  value  of  the  maximum  error  which  occurs  at  k  =  KkQ  can  be 
expressed  in  terms  of  as  follows: 
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The  corresponding  expression  for  the  error  in  the  function  kp(s)  due 

to  the  total  error  in  the  function  Fc(k)  shown  in  the  preceding  sketch 
i  S  ' 

“  26^^^  kl^8  )error  =  ~  _  ^  cOs  ^(^C)3)  ”  s^ko  Ci(iskos)  +  Si(iikQs)  - 


Si(kQs)  +  sin(Kkos) 


(B5) 


A  plot  of  equation  (B5)  as  a  function  of  IcqS  is  shown  in  figure  8 
for  several  values  of  the  parameter  k.  Thus,  for  a  maximum  error  emax 
in  the  Fc(k)  function  at  any  point  kQ  <  Kk q  <  «>,  an  estimate  of  the 

corresponding  error  in  the  function  kpCs)  can  he  obtained  from 
figure  8. 
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Figure  1.-  Comparison  of  indicial  lift  functions  due  to  a  sudden  change 
in  sinking  velocity  or  pitching  displacement  for  several  subsonic 


Asymptotic  values 
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Figure  3--  Comparison  of  indicial  lift  functions  due  to  a  sudden  change 
in  pitching  velocity  for  an  airfoil  rotating  about  its  leading  edge 
for  several  subsonic  Mach  numbers. 
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Figure  7*-  Comparison  of  the  indicial  lift  functions  due  to  penetration 
of  a  sharp-edge  gust  at  several  subsonic  Mach  numbers. 
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(a)  Sinking  velocity. 


(b)  Pitching  velocity. 

Figure  10.-  Comparison  of  in-phase  oscillatory  moment  coefficients  u* 
in  present  analysis  with  coefficients  determined  from  reference  7. 
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